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ABSTRACT 

Existing  analytical  models  of  wave-induced  minor  constituent  fluctuations  result  from  linearized  perturbation 
expansions  of  a  rate  equation  governing  either  number  density  or  mixing  ratio,  whereas  many  numerical  models 
now  use  isentropic  parcel  advection  methods  to  simulate  these  effects.  Exact  relationships  between  the  two 
approaches  are  not  currently  clear  for  gravity  waves.  Here,  the  parcel  advection  method  is  formalized  and  applied 
to  derive  analytical  formulas  for  the  response  of  vertical  constituent  profiles  of  arbitrary  shape  to  adiabatic 
gravity  wave  displacements.  These  relations  are  compared  to  corresponding  formulas  from  standard  linearized 
perturbation  analyses.  Both  methods  accurately  model  perturbations  produced  by  nondissipating  hydrostatic 
gravity  waves  within  idealized  vertical  tracer  profiles.  Both  methods  can  also  model  wave-induced  perturbations 
of  minor  constituents  with  shorter  chemical  lifetimes.  This  is  demonstrated  by  using  the  parcel  method  to 
reproduce  previous  results  for  wave-induced  fluctuations  in  upper-stratospheric  ozone.  The  parcel-based  approach 
yields  more  accurate  models  of  nondissipating  hydrostatic  gravity  wave  effects  on  tracer  distibutions  with  sharp 
spatial  gradients.  The  authors  demonstrate  this  by  using  the  method  to  model  wave-induced  modulations  of 
sporadic  sodium  layers  in  the  mesosphere  and  of  ozone  laminae  in  the  lower  stratosphere.  Conversely,  the  parcel 
method  does  not  accurately  model  tracer  density  perturbations  produced  by  nonhydrostatic  waves,  or  situations 
in  which  the  photochemical  response  of  the  constituent  leads  to  significant  diabatic  feedback  on  the  perturbing 
wave. 


1.  Introduction 

High-resolution  observations  of  minor  (trace)  con¬ 
stituents  in  the  middle  and  upper  atmosphere  often  re¬ 
veal  small-scale  perturbations  superimposed  upon  the 
mean  background  distribution.  Such  features  can  be  pro¬ 
duced  by  any  number  and  combination  of  dynamical, 
radiative,  and  chemical  processes.  For  trace  constituents 
with  long  lifetimes,  however,  dynamics  must  be  pre- 
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dominately  responsible  for  any  observed  variability  of 
that  constituent’s  local  number  density  over  shorter 
timescales. 

One  important  source  of  such  variability  is  waves. 
For  instance,  vertical  soundings  at  middle  and  high  lat¬ 
itudes  often  reveal  large  enhancements  or  depletions  of 
lower-stratospheric  ozone  (a  tracer  at  these  heights)  con¬ 
fined  to  narrow  (—1-2  km)  vertical  layers  (e.g.,  Reid 
and  Vaughan  1991;  Bird  et  al.  1997).  These  structures, 
known  assortedly  as  filaments  or  laminae,  result  from 
complex  synoptic-scale  flow  patterns  produced  by 
breaking  planetary  Rossby  waves,  which  yield  air  par¬ 
cels  at  different  stratospheric  levels  with  vastly  different 
histories  and,  hence,  constituent  densities  (e.g.,  Newman 
and  Schoeberl  1995;  Orsolini  et  al.  1997). 

Oblique  oscillatory  displacements  of  air  parcels  from 
their  equilibrium  positions  produce  accompanying  den¬ 
sity  and  temperature  changes  that  can  also  perturb  con¬ 
stituent  profiles.  Oscillations  due  to  various  types  of 
planetary  waves  have  been  found  to  induce  related  per¬ 
turbations  in  ozone  (e.g.,  Newman  and  Randel  1988; 
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Prata  1990;  Hess  1990;  Randel  and  Gille  1991;  Randel 
1993;  Stanford  and  Ziemke  1993;  Wirth  1993;  Ziemke 
and  Stanford  1994;  Engelen  1996)  and  other  trace  con¬ 
stituents  (e.g.,  Randel  1990;  Salby  et  al.  1990;  Lim- 
pasuvan  and  Leovy  1995).  On  smaller  scales,  internal 
gravity  waves  can  also  give  rise  to  fluctuations  in  con¬ 
stituent  profiles  (e.g.,  Reber  et  al.  1975;  Gardner  and 
Shelton  1985;  Hedin  and  Mayr  1987;  Sugiyama  1988; 
Hoegy  et  al.  1990;  Wilson  et  al.  1991;  Alexander  and 
Pfister  1995;  Langford  et  al.  1996;  Bian  et  al.  1996; 
Bacmeister  et  al.  1997). 

In  fact,  since  observational  data  on  middle  and  upper 
atmosphere  constituents  are  often  better  and  more  abun¬ 
dant  than  wind  and  temperature  data,  analyses  of  fluc¬ 
tuations  in  long-term  measurements  of  constituent  den¬ 
sities  have  been  used  to  infer  climatologies  of  both  plan¬ 
etary  waves  and  gravity  waves  (e.g.,  Randel  and  Gille 
1991;  Senft  and  Gardner  1991;  Randel  1993;  Ziemke 
and  Stanford  1994;  Collins  et  al.  1996).  However,  since 
the  response  of  constituents  to  waves  is  usually  not 
straightforward,  wave  properties  must  be  inferred  from 
such  observations  by  using  a  wave-tracer  interaction 
model  to  convert  the  measured  constituent  response  into 
a  wave-related  quantity.  Such  conversion  methods  have 
been  derived  to  date  by  applying  standard  linearized 
perturbation  methods  of  gravity  wave  theory  to  rate 
equations  governing  constituent  densities  or  mixing  ra¬ 
tios.  Such  derivations  yield  analytical  formulas  for 
translating  conventional  gravity  wave  oscillations  into 
a  corresponding  response  in  the  constituent  profile  (e.g., 
Dudis  and  Reber  1976;  Chiu  and  Ching  1978;  Gardner 
and  Shelton  1985).  Such  relations  can  then  be  used  to 
convert  observed  constituent  fluctuations  into  wave-re¬ 
lated  quantities  (e.g.,  Gardner  and  Voelz  1985;  Randel 
1990;  Senft  and  Gardner  1991). 

Recently,  detailed  numerical  models  have  provided 
new  insights  into  the  way  minor  constituents  respond 
to  planetary  wave  and  gravity  wave  advection  (e.g., 
Fritts  et  al.  1993,  1997;  Waugh  et  al.  1994;  Roble  and 
Shepherd  1997).  Rich  complex  structures  in  minor  spe¬ 
cies  distributions  are  simulated,  which  often  resemble 
observations  (e.g.,  Fritts  et  al.  1993,  1997)  but  are  not 
predicted  by  standard  formulas  that  convert  between 
observed  tracer  variability  and  a  plane  monochromatic 
wave.  Due  to  their  complexity,  however,  these  models 
cannot  as  yet  be  used  for  operational  extraction  of  wave 
properties  from  observed  constituent  fluctuations  in  the 
atmosphere.  On  the  other  hand,  simplified  versions  of 
these  models,  which  retained  the  minimum  terms  needed 
to  simulate  the  basic  effects,  might  lead  to  improved 
analytical  models  of  wave-induced  perturbations  of  mi¬ 
nor  constituent  profiles. 

We  investigate  this  possibility  here  by  considering 
one  class  of  these  numerical  models,  the  so-called  tra¬ 
jectory  or  parcel  advection  models.  They  have  proved 
successful  in  simulating  finescale  tracer  structures  in  the 
stratosphere,  often  using  wind  and  temperature  data  with 
much  coarser  spatiotemporal  resolution  (e.g.,  Newman 


and  Schoeberl  1995;  Orsolini  et  al.  1997).  The  models 
work  well  because  their  assumptions  of  adiabatic  motion 
and  long  chemical  lifetimes  are  well  satisified  for  many 
constituents  in  the  lower  stratosphere,  where  diabatic 
heating  and  turbulent  diffusivities  are  both  small.  How¬ 
ever,  constituents  with  shorter  chemical  lifetimes  can 
also  be  simulated  by  these  models  if  a  suitable  parcel- 
based  chemistry  scheme  is  used  (e.g.,  Lutman  et  al. 
1997).  Even  in  the  troposphere,  where  diffusion  and 
diabatic  heating  can  be  large,  parcel  advection  models 
that  parameterize  these  effects  have  also  proved  to  be 
useful  tools  (e.g.,  Kao  et  al.  1995;  Rind  and  Lerner 
1996). 

Given  that  parcel  advection  models  often  perform 
well  even  after  many  days  of  integration,  the  same  mod¬ 
els  should  also  perform  well  if  the  synoptic-scale  ad¬ 
vection  patterns  were  supplemented  with  wave-induced 
advection  effects  since  their  timescales  (periods)  are  of 
the  order  of  days  for  planetary  waves  and  less  than  a 
day  for  gravity  waves.  Furthermore,  since  many  gravity 
waves  and  planetary  waves  propagate  conservatively 
through  the  stratosphere,  an  assumption  of  adiabatic 
wave-induced  advection  should  also  be  well  satisfied. 
Indeed,  limited  tests  of  gravity  wave  effects  in  such 
models  have  used  this  sort  of  approach  (Pierce  et  al. 
1994;  Jensen  and  Thomas  1994;  Carslaw  et  al.  1998). 

Thus,  we  formalize  the  basic  physical  concepts  un¬ 
derpinning  parcel  advection  models  in  section  2,  assess 
the  utility  of  the  method  in  gravity  wave  studies,  then 
apply  it  in  section  3  to  derive  analytical  responses  of 
tracers  to  wave-induced  displacements.  We  compare  our 
results  throughout  with  corresponding  results  using  the 
traditional  approach  of  expanding  state  parameters  into 
the  sum  of  a  mean  and  perturbation  term,  then  isolating 
and  linearizing  the  perturbation  terms  in  the  governing 
equations.  Results  of  these  analyses  and  intercompari¬ 
sons  are  summarized  in  section  4. 

2.  Air  parcel  approach 

This  approach  is  based  on  the  classic  meteorological 
problem  of  vertically  displacing  a  small  dry  parcel  of 
air  from  its  equilibrium  height,  where  it  is  “well  mixed” 
and  in  thermal  equilibrium  with  the  surrounding  at¬ 
mosphere.  We  assume  that  the  displacement  is  small 
and  that  the  parcel  responds  adiabatically  while  always 
maintaining  a  pressure  equal  to  that  of  the  surrounding 
atmosphere,  following  standard  conventions  (see,  e.g., 
section  2.5.1  of  Wallace  and  Hobbs  1977;  section  9.1 
of  Iribarne  and  Godson  1981). 

a.  Conserved  quantities  in  the  parcel 

1)  Tracer  mixing  ratio 

We  consider  a  constituent  of  number  density  n  within 
an  air  parcel  of  total  number  density  nM.  For  a  trace 
constituent  n/nM  <C  1,  any  changes  in  n  produce  neg- 
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ligible  changes  in  nM,  whereupon  the  following  conti¬ 
nuity  equations  apply: 

+  nMV  •  U  =  0  (la) 

at 

dn  _ 

—  +  rcV  •  U  =  R,  (lb) 


where  U  is  the  local  velocity  vector,  d/dt  =  d/dt  +  U  •  V 
is  a  time  derivative  in  the  Lagrangian  (parcel  following) 
frame  (d/dt  is  the  ground-based  Eulerian  time  deriva¬ 
tive),  and  R  is  the  net  production/loss  rate  for  the  trace 
constituent’s  number  density  n.  Eliminating  V  •  U  be¬ 
tween  (la)  and  (lb)  yields  (Lindzen  and  Goody  1965) 


d  (  n\  _  dq  _  R 
dt\nM)  dt  nM' 


(2) 


where  q  is  the  mixing  ratio  of  the  trace  constituent. 

For  passive  trace  constituents,  R  ~  0,  whereupon  (2) 
yields 


dq 

dt 


=  0, 


(3a) 


or 


n 

q  =  —  =  const  (3b) 

«M 

within  the  parcel.  Thus,  we  get  the  well-known  result 
that  passive  trace  constituents  are  also  tracers  since  the 
movement  of  a  region  of  unique  constituent  mixing  ratio 
traces  the  Lagrangian  motion  of  air  parcels. 


2)  Potential  temperature 


Air  parcels  advected  by  an  inviscid  adiabatic  flow 
also  conserve  potential  temperature  ©  so  that 


(4) 


For  a  horizontally  homogeneous  atmosphere,  the  back¬ 
ground  potential  temperature  ©  is  given  by  Poisson’s 
relation 


/— /  \\  (y-i)/y 

©(z)  =  =  T(z,  z0),  (5) 

where  z  is  height,  y  is  the  ratio  of  specific  heats,  T(z) 
is  the  background  temperature,  and  ~p(z)  is  the  back¬ 
ground  pressure  (e.g..  Turner  1973).  From  (5),  © (z)  is 
also  equal  to  the  temperature  T(z,  z0)  that  an  air  parcel 
attains  on  being  transported  adiabatically  a  distance  Sz 
from  its  rest  height  z  to  a  new  height  z0  =  z  +  Sz. 

3)  Potential  density 

We  can  also  work  in  terms  of  potential  density,  D(z), 
which  is  defined  for  a  horizontally  homogeneous  back¬ 
ground  atmosphere  as 


D(z)  =  p(z)^^j  =  p(z,  z0)  (6) 

(e.g.  Turner  1973),  where  JHz.)  =  MnM(z)  is  the  back¬ 
ground  total  air  density  at  height  z,  nM(z)  is  the  back¬ 
ground  total  number  density,  and  M  is  the  mean  mass 
of  an  air  molecule,  assumed  constant  (no  diffusive  sep¬ 
aration).  From  (6),  D(z )  also  equals  the  density  p(z,  z0) 
that  an  air  parcel  attains  on  being  transported  adiabat¬ 
ically  from  a  pressure,  p(z)  and  density  JHz),  at  its  rest 
height  z  to  a  reference  pressure  p(z0)  at  the  displaced 
height  z0  =  z  +  Sz.  It  then  follows  (e.g.,  Wallace  and 
Hobbs  1977)  that 


^  =  0 
dt 


(7) 


for  parcels  on  an  inviscid  adiabatic  flow,  and  thus  the 
potential  density  of  any  parcel  within  such  a  flow  is  also 
conserved. 


4)  Potential  tracer  density 

Analogously,  we  can  also  define  a  potential  density, 
A,  for  a  given  tracer.  If  the  constituent  in  question  be¬ 
haves  as  a  tracer  at  all  heights  and  has  negligible  hor¬ 
izontal  gradients  in  its  distribution,  then,  as  in  (6),  we 
can  define  a  background  potential  tracer  density  profile 
given  by 


A(z)  =  Mlcv(z )  =  Mtcn(z)(^y^]j  =  Mtcn(z,  z0), 

(8) 

where  Mtc  is  the  molecular  mass  of  the  trace  constituent 
and  ~v(z)  is  the  background  potential  number  density  of 
the  tracer.  From  (8),  we  see  that  A(z)  also  equals  the 
mass  density  of  tracer  Mlcn(z,  z0)  that  an  air  parcel  con¬ 
tains  after  being  transported  adiabatically  from  a  pres¬ 
sure  p(z)  and  tracer  density  Mlcn(z )  at  its  rest  height  z 
to  a  reference  pressure  p(z0)  at  the  new  height  z„  =  z 
+  Sz.  It  then  follows  that 


d  A  dv 
dt  dt 


(9) 


for  parcels  on  an  inviscid  adiabatic  flow.  Since  number 
densities  are  more  commonly  considered  in  tracer  stud¬ 
ies,  we  shall  work  in  terms  of  v  rather  than  A  in  sub¬ 
sequent  analysis. 


b.  Perturbations  due  to  adiabatic  vertical 
displacements 

Next  we  consider  local  responses  of  each  conserved 
term  to  a  small  vertical  adiabatic  displacement  £  of  an 
air  parcel  from  its  equilibrium  position  z,  to  a  new  height 

z2  =  z,  +  i,  (10) 

given  the  standard  background  pressure  profile 
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(a)  Potential  Temperature 


(c)  Potential  Density 


P 


Fig.  1.  Schematic  diagram  of  response  (light  dotted  curve)  of  (a)  potential  temperature  ©,  (b)  absolute  temperature  T,  (c)  potential  density 
D,  and  (d)  absolute  density  p  in  an  air  parcel  (shown  as  a  cube)  as  it  is  adiabatically  advected  a  small  vertical  distance  £  from  its  equilibrium 
rest  position  zx,  where  it  is  in  thermal  equilibrium  with  the  background  atmosphere.  Background  profiles  of  each  quantity  are  plotted  with 
solid  curves. 


p  fe) 


p(z  i) 


“  P(Zi) 


(11) 


where  H(zt)  =  RT{zi)lg  is  the  pressure  scale  height  at 
height  Zu  and  the  latter  expression  in  (1 1)  is  sufficiently 
accurate  for  small  £  given  a  background  temperature 
profile  T(z)  that  varies  slowly  with  height. 

The  formal  parcel-based  analysis  of  this  problem  is 
given  in  appendix  A,  where  departures  (perturbations) 
of  parcel-based  quantities  from  surrounding  background 
values  are  derived.  The  basic  procedure  for  computing 
these  perturbations  is  depicted  schematically  in  Fig.  1 
for  both  temperature  and  density.  In  each  panel  in  Fig. 
1,  an  air  parcel  has  been  advected  from  its  equilibrium 
height  Zi  to  a  new  height,  z2,  and  in  each  case  the  adi¬ 
abatic  response  is  given  by  moving  the  parcel  along  the 
light  dotted  “response  curve.”  The  top  panels  of  Fig. 
1  show  constructions  for  the  potential  temperature  and 


potential  density:  since  both  remain  constant  within  the 
air  parcel,  their  response  curves  are  vertical.  Corre¬ 
sponding  constructions  in  terms  of  the  absolute  tem¬ 
perature  and  density  are  shown  in  the  lower  panels. 
Their  response  curves  are  given  by  T(z2,  z2)  and  p(z, , 
z2)  in  (5)-(6):  hereafter,  it  proves  expedient  to  use  the 
abbreviated  notation  T(z2)  and  p(z2)  for  these  expres¬ 
sions  since  z,  stays  constant  and  the  functional  depen¬ 
dence  on  z,  is  implicit  from  (10).  Background  profiles 
are  shown  with  thick  solid  curves  on  each  plot.  Per¬ 
turbations  are  given  by  the  difference  between  the  adi¬ 
abatically  varying  value  within  the  parcel  and  the  local 
background  value  at  height  z2,  as  shown  in  Fig.  1.  Using 
the  abbreviated  notation  introduced  above,  this  yields 
perturbation  terms  ©'(Z2X  D'(z2),  T'(z2),  and  p'{z2)- 
Simplified  equations  for  these  perturbation  terms  fol¬ 
low  when  the  displacement  £  is  sufficiently  small  that 
variations  in  both  the  background  profiles  and  the  ad¬ 
iabatic  response  curves  between  z.t  and  z.2  are  approx¬ 
imately  linear  in  Fig.  1.  This  holds  so  long  as  any  cur¬ 
vature  of  these  profiles  occurs  over  vertical  scales  L  S> 
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£  Since  the  formal  parcel  relations  in  appendix  A  are 
height-varying  exponentials,  L  can  be  conveniently 
characterized  by  the  scale  height  Hx  of  the  exponential 
relation  for  quantity  X,  so  that  I  zf///Yl  <C  1.  When  this 
condition  is  satisfied,  the  exponential  parcel  relations 
can  be  accurately  approximated  by  retaining  only  the 
first  two  terms  from  their  MacLaurin  series  expansions. 
For  temperatures  and  densities,  these  truncated  Mac¬ 
Laurin  series  (TMS)  relations  yield  the  instantaneous 
simplified  perturbation  relations 


p'fe)  =  £fe)  =  /B'feA  =  (r(z2)\ 

pfe)  D(z2)  \©fe)/  \rfe)J 


=  C  N2(Zj) 

Hd(z i)  g 

where 


(12) 


Hd(z  i) 


mzjY1 

dz  ) 


g 

N2(z ,) 


(13) 


is  the  potential  density  scale  height  (see  appendix  A). 
Equation  (12)  holds  for  l£///Dl  <C  1. 


c.  Range  of  applicability  to  gravity  waves 

Next  we  consider  a  sinusoidal  form  for  the  vertical 
parcel  displacement  from  z, : 

£  =  C(x,  y,  zt,  t ) 

=  l  sin(kx  +  ly  +  mzt  —  art  +  cp),  (14) 

where  %  is  the  peak  vertical  displacement  amplitude,  cp 
is  a  (constant)  phase  offset,  ( k ,  l,  in)  is  the  wavenumber 
vector,  and  co  is  the  oscillation  frequency  of  the  parcel 
about  its  equilibrium  height  z,  .  Considering  potential 
temperature  only  for  the  moment,  (14)  and  (12)  yield 


The  relations  (14)  and  (18)  have  the  same  form  as 
the  simplified  polarization  relations  that  govern  a  hy¬ 
drostatic  gravity  wave  of  intrinsic  frequency  a>.  This 
occurs  because,  in  the  linear  hydrostatic  limit,  a  gravity 
wave  advects  air  parcels  adiabatically  (e.g.,  Fritts  and 
Rastogi  1985),  has  a  small  displacement  amplitude 
l  (I £/Hd\  <C  1),  and  produces  a  negligible  pressure  per¬ 
turbation  within  displaced  parcels  (e.g.,  Gossard  and 
Hooke  1975). 

To  assess  the  degree  and  range  of  this  correspondence 
more  quantitatively,  a  standard  non-Boussinesq  linear¬ 
ized  perturbation  derivation  of  the  acoustic-gravity 
wave  equations  is  outlined  in  appendix  B.  There  it  is 
shown  that  the  parcel-derived  polarization  relations  (18) 
are  an  acceptable  approximation  to  the  full  acoustic- 
gravity  wave  polarization  relations  when  the  normalized 
quantities  I a(m,  co) I  and  \b(m,  w)l  are  <Cl  [see  Eqs.  (BIO) 
and  (B 11)].  When  this  limit  is  not  satisfied,  wave-in¬ 
duced  pressure  perturbations  become  nonnegligible  [see 
Eq.  (B 1 3)],  so  the  parcel  approach  (which  assumes  zero 
pressure  perturbation)  becomes  inaccurate.  Sample  cal¬ 
culations  in  appendix  B  show  that  the  differences  be¬ 
tween  the  full  acoustic-gravity  wave  polarization  re¬ 
lations  and  the  parcel-based  approximations  (18)  are 
small  for  \z  =  2tt/\iii\  &  20  km  (see  Figs.  Bib  and  Blc). 
Thus,  we  use  this  as  an  approximate  upper- wavelength 
bound  on  the  validity  of  the  hydrostatic  gravity  wave 
approximation  and  of  subsequent  parcel-based  deriva¬ 
tions  of  gravity  wave-induced  minor  constituent  vari¬ 
ability. 

3.  Gravity  wave  perturbations  of  minor 

constituents 


Q'fo  y,  z2,  t) 
©  (z2) 


-NHzJ 

g 


£'(x,  y,  zu  t). 


(15) 


Since  l  is  small  (I  £/HD\  <  1),  replacing  the  variable 
term  ©(z2)  with  the  constant  value  @(z,)  introduces 
negligible  errors  in  (15).  Then,  to  a  good  approximation, 
we  can  remove  the  explicit  transport  of  parcels  to  z2  so 
that  we  can  reexpress  the  potential  temperature  pertur¬ 
bation  (15)  as  a  sinusoidal  expression  of  the  form  (14) 
at  a  fixed  height  z, : 

©'(a,  y,  Zi,t)  =  @(m)  sin  (lac  +  ly  +  mzt  —  cot  +  cp),  (16) 


where 


©(m) 


~N2(zt) 

g 


©(*,)£ 


d@fei)i 
dz  ^ 


(17) 


We  now  apply  the  parcel-based  approach  to  model 
gravity-wave-induced  perturbations  of  minor  constitu¬ 
ent  profiles.  We  begin  by  considering  some  simple  ex¬ 
amples  and  compare  the  parcel-based  derivations  (ap¬ 
pendix  A)  with  corresponding  equations  derived  by  a 
conventional  linearized  perturbation  analysis  (appendix 
B). 

a.  Tracer  with  linear  q(z)  profile 

The  parcel-based  approach  to  this  problem  is  depicted 
schematically  in  Fig.  2a  and  yields  the  perturbation  for¬ 
mula  (A  19).  Since  the  background  profile  is  linear,  then 
q(z2)  =  q(Zi)  +  ( dq/dz)£ ',  whereupon  (A19)  and  the 
wave  solution  (14)  yield 


Transport  of  parcels  to  z2  is  now  implicit  within  (16) 
and  (17).  Similar  relations  to  (17)  follow  from  (12)  and 
(14)  for  relative  perturbations  of  potential  density,  den¬ 
sity,  and  temperature,  with  peak  amplitudes  interrelated 
as 


q'(x,  y,  z2,  t)  =  ~  y,  Zj,  t ).  (19) 

As  for  (15),  if  l  is  small,  then  the  explicit  transport  of 
parcels  to  z2  can  be  removed  in  order  that  (19)  can  be 
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(a)  Tracer  Mixing  Ratio 


q 


(b)  Tracer  Number  Density 


(c)  Potential  Tracer  Number  Density 


Fig.  2.  As  in  Fig.  1  but  for  (a)  tracer  mixing  ratio  q,  (b)  tracer 
number  density  n,  and  (c)  potential  tracer  number  density  v  of  the 
parcel. 


reexpressed  as  a  wave  solution  of  the  form  (16),  with 
peak  amplitudes  related  as 


(20) 

A  standard  linearized  perturbation  analysis  of  (3a) 
using  the  methodology  in  appendix  B  yields  the  lin¬ 
earized  perturbation  mixing  ratio  equation  (e.g.,  Randel 
1990) 


—  uoq  + 


(21) 


which,  since  w  =  —  also  yields  (20)  at  a  given 
height  z,  .  Thus  the  two  approaches  give  identical  results 
in  this  case. 


b.  Tracer  with  a  linear  n(z)  profile 

The  parcel-based  approach  to  this  problem  is  depicted 
schematically  in  terms  of  tracer  number  density  in  Fig. 
2b  and  in  terms  of  potential  tracer  number  density  in 
Fig.  2c.  The  mathematical  derivations  are  set  out  in 
appendix  A.  Because  n(z)  is  linear  and  v{z)  is  approx¬ 
imately  linear  over  the  interval  £  then  TMS  expansions 
of  the  parcel  relations  for  tracer  perturbations  [e.g., 
(A25b)]  are  accurate,  leading  to  the  simplified  TMS 
result 


v'{x,  y,  z2,  t)  _  n'(x,  y,  z2 ,  0 
v(z2)  n{z2 ) 


C(x,  y,  zu  t) 


HV{Z\) 


(22) 


[see  (A26)],  where 


H£zd 


(  1  +  1  anfafr1 

lytf(Zi)  «(Zi)  dz  ) 


(23) 


is  the  “potential  tracer  density  scale  height”  by  analogy 
to  the  potential  total  density  scale  height  (13)  [see  (A7)]. 
Since  ii(z)  is  linear,  the  accuracy  of  the  TMS  relation 
(22)  is  controlled  only  by  the  curvature  of  the  expo¬ 
nential  response  curve  (A21)  shown  in  Fig.  2b  and  so 
is  accurate  for  \%/yH\  <C  1. 

As  for  (15)  and  (T9),  if  l  is  small,  then  the  explicit 
transport  of  parcels  to  z2  can  be  removed  from  (22)  so 
that  this  equation  can  be  reexpressed  as  a  wave  solution 
of  the  form  (16)  at  the  equilibrium  height  z1,  with  peak 
amplitudes  related  as 


n(Zi)  _  11  dn 

«(Zi)  yH{zd  n(Zi )  dz _ 

where  dii/dz  is  a  constant  here. 

Equation  (24)  also  follows  from  a  linearized  pertur¬ 
bation  analysis  using  the  nonhydrostatic  isothermal 
acoustic-gravity  wave  equations  of  Hines  (1960)  (Dudis 
and  Reber  1976;  Chiu  and  Ching  1978;  Weinstock  1978; 
Gardner  and  Shelton  1985;  Sugiyama  1988).  A  similar 
analysis  in  appendix  B  shows  that  the  full  expression 
for  nonhydrostatic  acoustic-gravity  waves  is  [see  (B20)] 


l  (24) 


n(z,)  _  I  1  1  dn 

«(Zi)  1  ytf(Zi)  +  n(Zi)  dz 


N2(zfi 

8 


[a(m,  co)  +  ib{m,  w)] 


l 


(25) 


Therefore  (24)  is  valid  only  when  I a(m,  <w)l  and  I b(m, 
co)  I  are  negligible,  which,  as  discussed  in  section  lc  and 
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Nonlinear  n(z)  Profile 


Fig.  3.  As  in  Fig.  2b  but  illustrating  the  complications  introduced 
by  advecting  a  parcel  across  a  nonlinear  n(z )  profile. 


appendix  B,  holds  only  for  hydrostatic  gravity  waves 
of  A_  S  20  km.  Thus,  (24)  is  a  hydrostatic  result  (see 
also  Dudis  and  Reber  1976),  despite  the  fact  that  non¬ 
hydrostatic  wave  equations  are  needed  to  arrive  at  this 
result  in  some  linearized  perturbation  analyses  (e.g., 
Gardner  and  Shelton  1985).  Note  that  no  such  caveats 
applied  to  the  mixing  ratio  relation  (20) — it  holds  for 
both  hydrostatic  and  nonhydrostatic  waves.  Therefore, 
mixing  ratio  is  the  more  straightforward  quantity  to  con¬ 
sider  when  converting  between  minor  constituent  per¬ 
turbations  and  displacement  perturbations  due  to  non¬ 
hydrostatic  gravity  waves  or  “tall”  (fast)  Kelvin  waves. 


c.  Tracer  with  a  weakly  nonlinear  n(z)  profile 


In  the  parcel-based  approach,  the  assessment  is 
straightforward.  For  an  exponential  n{z)  profile,  for  in¬ 
stance,  (22)  and  hence  (24)  are  accurate  for  I  %/HJ  <C 
1,  from  (A6a)  and  (A6b).  For  arbitrary  profile  shapes, 
we  simply  avoid  TMS  approximations  and  use  the  gen¬ 
eral  parcel  solution  (A21),  which,  using  (14),  yields  the 
normalized  perturbation  relation 


n'(x,  y,  z2,  t ) 
n  (z2) 


«(m)  (~C(x,  y,  Zi,  t)\ 

^fe)eXPl  T/Zfc)  j 


-  1, 


(26) 


from  (A23b),  where  z2  ==  Z,  +  C(x<  y>  Z<  ,  0-  Note  that 
(26)  yields  a  nonsinusoidal  solution.  In  terms  of  poten¬ 
tial  tracer  density  v,  this  is  equivalent  to  using  (A6a) 
or  (A8)  instead  of  the  approximate  relation  (A6b).  Direct 
comparison  of  results  using  (24)  and  (26)  then  indicate 
whether  the  former  equation  is  sufficiently  accurate  for 
modeling  wave  perturbations  of  a  given  n(z)  profile. 

Such  assessments  are  less  straightforward  within  the 
framework  of  a  linearized  perturbation  analysis  (e.g., 
Chiu  and  Ching  1978).  Gardner  and  Shelton  (1985)  de¬ 
veloped  an  iterative  method  for  deriving  higher-order 
corrections  to  (24)  due  to  the  curvature  of  a  given  n{z) 
profile,  using  perturbation  expansions  based  on  an  as¬ 
sumed  form  for  the  solution.  Their  most  general  solution 
took  the  form 


nix,  y,  z,  t) 


z  ~  yH(z)  In  1  + 


A'(x,  y,  z,  t) 

Y  -  1 


1  + 


A'(x,  y,  z,  t) 
y  -  1 


(27) 


Of  course  tracer  density  profiles  cannot  be  infinitely 
linear,  but  instead  usually  approximate  some  type  of 
Chapman  curve.  If  the  curvature  of  the  profile  is  gradual 
with  height,  then,  for  hydrostatic  waves,  (24)  is  often 
used  to  model  wave-induced  perturbations,  with  dnizfil 
dz  reevaluated  at  each  new  height  in  (24),  as  in  (23) 
(e.g.,  Chiu  and  Ching  1978;  Gardner  and  Voelz  1985; 
Senft  and  Gardner  1991).  For  this  to  be  accurate,  the 
profile  needs  to  be  approximately  linear  only  over  the 
limited  range  Zi  ±  £,  as  shown  in  Fig.  2b.  However, 
this  approximation  may  not  always  be  sufficiently  ac¬ 
curate  for  certain  profile  shapes,  as  shown  schematically 
in  Fig  3.  Thus,  tests  of  its  adequacy  for  a  given  profile 
are  required. 


where  A'(x,  y,  z,  t)  is  a  normalized  gravity  wave-induced 
perturbation,  which,  for  small-amplitude  waves  in  an 
isothermal  atmosphere  (Hines  1960),  is  approximately 
equal  to  the  relative  density  perturbation  of  the  wave 
[see  Eq.  (38)  of  Gardner  and  Shelton  1985].  It  then 
follows  that  (27)  corresponds  approximately  to  the  par¬ 
cel  solution  (A21)  derived  in  appendix  A  since  1  + 
A'{x,  y,  z,  t)/(y  —  1)  ~  1  +  £'{x,  y,  z,  t)/yH(z),  which 
in  turn  is  a  TMS  expansion  of  exp[£'(x,  y,  z,  t)/yH(z)\. 

Gardner  and  Shelton  converted  their  solutions  into 
perturbation  formulas  at  a  set  height  z,  using  MacLaurin 
series  expansions.  They  derived  the  following  extension 
of  (24)  in  which  the  effects  of  both  first-  and  second- 
order  vertical  derivatives  of  n{z)  were  retained: 


n'(x,  y,  zu  t) 

1  + 
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(28) 
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where  p'(x,  y,  z, ,  t)  is  the  gravity  wave  density  oscil¬ 
lation,  which  in  this  case  is  related  to  l  through  (14) 
and  (18)  using  isothermal  relations  (dH(zx)ldz  =  0)  to 
evaluate  N2(zi)  in  (13).  Like  (26),  Eq.  (28)  is  also  a 
nonsinusoidal  solution. 

In  Fig.  4,  we  compare  (28)  with  the  linear  profile 
solution  (24)  and  the  parcel  solution  (26).  Since  sodium 
acts  as  a  passive  tracer  to  gravity  waves  above  ~85  km 
(e.g.,  Hickey  and  Plane  1995),  we  use  an  n(z )  profile 
based  on  a  Gaussian  model  of  the  mesospheric  sodium 
layer,  similar  to  the  one  considered  by  Gardner  and  Shel¬ 
ton  (1985).  The  profile  is  plotted  as  the  thick  light  curve 
in  Figs.  4b  and  4e.  We  perturb  this  profile  with  a  hy¬ 
drostatic  gravity  wave  (14)  of  vertical  wavelength  A_  = 
10  km,  l  =  800  m  (both  representative  mesospheric 
values),  and  constant  phase  <p0 .  Sample  waves  are  shown 
in  Figs.  4a  and  4d  as  a  function  of  height  z  at  a  given 
location  (x0,  y0)  at  two  different  times  t.  The  corre¬ 
sponding  perturbations  to  the  background  tracer  profile 
are  plotted  in  Figs.  4b  and  4e,  and  the  normalized  per¬ 
turbations  are  isolated  and  plotted  in  Figs.  4c  and  4f  for 
the  parcel  solution  (26)  and  the  perturbation  formulas 
(24)  and  (28).  The  parcel  solution  (26)  is  also  replotted 
with  a  dotted  curve  in  Figs.  4c  and  4f  at  the  undisturbed 
height  z  =  Zi,  instead  of  at  the  perturbed  height  z  =  z2 
=  z,  +  £(.x0,  y0,  z, ,  t )  (solid  curve),  to  help  separate 
the  deviations  that  arise  due  to  different  physical  lo¬ 
cations  for  each  parcel  (i.e.,  differences  in  z2  and  Z\ ) 
and  those  due  to  incomplete  characterization  of  the  ef¬ 
fects  of  mean  profile  curvature. 

We  see  in  Figs.  4c  and  4f  that  the  simple  relation  (24) 
is  reasonably  accurate  here  but  that  the  higher-order 
correction  (28)  brings  the  solution  closer  to  the  parcel 
solution  (26).  The  results  in  Fig.  4  also  show  many  well- 
known  features  of  wave-induced  fluctuations  of  tracer 
number  density  profiles;  for  example,  that  the  wave  re¬ 
sponse  is  always  larger  on  the  underside  of  the  layer 
than  on  the  topside  (e.g.,  Chiu  and  Ching  1978;  Wein- 
stock  1978).  This  feature  is  easily  appreciated  in  a  parcel 
analysis  by  noting  that  all  points  on  the  unperturbed 
Gaussian  profile  in  Figs.  4b  and  4e  follow  the  “adiabatic 
rarefaction”  response  curve  (A21)  when  vertically  ad- 
vected,  as  shown  in  Fig.  2b.  This  clearly  yields  greater 
deviations  from  the  background  profile  on  the  bottom 
side  of  the  layer  than  on  the  top  side.  Also  evident  in 
Figs.  4c  and  4f  is  the  well-known  half-wavelength  phase 
flip  across  the  layer  peak  (Gardner  and  Shelton  1985). 

d.  Tracer  with  a  strongly  nonlinear  n(z)  profile 

Tracer  layers  in  the  middle  atmosphere  often  deviate 
significantly  from  the  smoothly  curved  profiles  consid¬ 


ered  in  Fig.  4.  For  example,  the  mesospheric  sodium 
layer  sometimes  develops  a  large  narrow  spike  of  en¬ 
hanced  Na  density,  known  as  a  sudden  (or  sporadic) 
sodium  layer  (SSL)  (Clemesha  et  al.  1980;  Kwon  et  al. 
1988;  Nagasawa  and  Abo  1995).  Similarly,  ozone  pro¬ 
files  in  the  lower  stratosphere  near  the  edge  of  the  polar 
vortex  often  contain  large  narrow  enhancements  or  de¬ 
pletions,  known  as  “laminae”  (e.g.,  Reid  and  Vaughan 
1991).  Neither  the  linearized  perturbation  relation  (24) 
nor  its  second-order  correction  (28)  are  sufficiently  ac¬ 
curate  for  modeling  wave-induced  effects  within  such 
profiles. 

However,  parcel-based  methods  can  be  used  since 
they  require  no  prior  assumptions  about  the  shape  of 
the  background  vertical  profile.  To  illustrate  this,  we 
simulate  an  SSL  by  adding  a  narrow  Gaussian  layer  to 
the  original  Gaussian  in  Fig.  4,  choosing  its  parameters 
to  resemble  the  “typical”  SSL  considered  by  Cox  et  al. 
(1993).  The  resultant  profile  is  the  thick  light  curve 
shown  in  each  panel  of  Fig.  5.  We  perturb  this  profile 
with  the  same  wave  as  in  Fig.  4.  Solid  curves  in  Fig. 
5  show  the  perturbed  profile  as  calculated  using  the 
parcel  solution  (A21)  and  the  wave  oscillation  (14),  with 
each  panel  showing  results  at  six  successive  equispaced 
time  values  1  spanning  a  full  wave  period  2  tt! oj.  We  see 
that  the  SSL  is  not  only  advected  up  and  down  by  the 
wave,  but  is  also  periodically  expanded  (Fig.  5b)  and 
compressed  (Fig.  5e). 

Detailed  features  of  the  response  are  more  easily  seen 
in  surface  renderings  of  successive  profiles,  shown  in 
Fig.  6.  Figure  6a  shows  the  periodic  oscillation  of  the 
SSL  peak  and  the  downward  movement  of  wave  phase 
fronts  above  and  below  the  SSL,  the  latter  associated 
with  an  approximately  linear  wave  response  within  the 
basic  Gaussian  profile  (Figs.  4c  and  4f).  The  SSL  is 
focused  upon  in  Fig.  6b,  which  shows  a  sinusoidal  ad- 
vection  of  its  peak  value  and  periodic  contraction  and 
expansion  of  its  width. 

Figure  7  shows  a  representative  lower-stratospheric 
ozone  profile  with  a  lamina  superimposed  (thick  light 
curve),  based  on  the  lidar  observations  of  Gibson-Wilde 
et  al.  (1997).  Like  sodium,  ozone  has  a  long  chemical 
lifetime  at  these  heights  and  therefore  acts  as  a  tracer 
to  gravity  wave  motions.  We  perturb  the  profile  with  a 
gravity  wave  of  £  =  200  m  and  Az  =  2.3  km,  choices 
based  on  the  gravity  wave  detected  and  characterized 
by  Gibson-Wilde  et  al.  (1997)  at  —  1 2—  1 8  km  during 
their  ozone  measurement.  Parcel  results  using  (A21)  for 
a  full  wave  cycle  are  shown  with  solid  curves  in  Fig. 
7.  Again,  we  see  not  only  vertical  advection  of  the  1am- 


— > 

Fig.  4.  Perturbations  of  a  Gaussian  number  density  profile  n(z )  [thick  light  curve  in  (b)  and  (e)]  are  plotted  in  (b)  and  (e),  as  produced 
by  the  gravity  wave  shown  above  it  in  (a)  and  (d),  respectively.  Corresponding  relative  perturbations  are  plotted  below  in  (c)  and  (f ).  Results 
using  parcel  advection  methods  (solid  and  dotted  curves)  are  compared  to  the  results  (24)  (dashed  curve)  and  (28)  (broken  curve)  based  on 
linearized  perturbation  solutions.  The  thick  light  curves  in  (c)  and  (f )  show  the  asymmetric  amplitude  envelope  of  the  parcel  solution,  given 
from  (A23a)  by  [n(z)!n{z  —  |)]  exp(  +  £/y//)  —  1.  An  isothermal  atmosphere  of  H  =  7  km  was  used. 


height  (km)  height  (km)  ^  height  (km) 
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(a)  Vertical  Displacement  Wave  C(xa,y0,z,,t)  (d)  Vertical  Displacement  Wave  £(x0,y0,Zj,t) 


Wave-Perturbed  Profiles  of  Tracer  Number  Density 


(e)  Wave-Perturbed  Profiles  of  Tracer  Number  Density 


Relative  Perturbations  of  Tracer  Number  Density 


(f)  Relative  Perturbations  of  Tracer  Number  Density 
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(a)  cot  =  0° 


(d)  cot  =  180" 


(b)  cot  =  60" 


(e)  cot  =  240" 


(c)  cot  =  120° 


(f)  cot  =  300° 


Fig.  5.  Sequence  of  six  profiles  (at  equispaced  times  t  spanning  one  wave  period)  of  a  wave-perturbed  Na  layer  (solid  curve)  as  in  Fig.  4, 
but  with  an  SSL  inserted  at  95  km  (based  on  Fig.  1  of  Cox  et  al.  1993).  The  thick  light  curve  shows  the  unperturbed  profile  in  each  case. 


ina,  but  compression  and  expansion  of  its  width  as  well, 
features  that  are  more  obvious  in  the  corresponding  sur¬ 
face  plot  in  Fig.  8.  Effects  such  as  these  may  account 
for  some  of  the  standard  deviation  in  mean  magnitudes 
and  widths  of  ozone  laminae  found  in  collated  ozone- 
sonde  data  (e.g..  Fig.  7  of  Reid  and  Vaughan  1991),  as 
well  as  variability  of  other  narrow  atmospheric  layers: 
for  example,  tropospheric  chemical  layers  (Newell  et  al. 
1996),  cloud  structures  (e.g.,  Carslaw  et  al.  1998),  polar 
mesosphere  summer  echoes  (e.g..  Fig.  1  of  Rtister  et  al. 
1996),  and  stratospheric  aerosols  (Shibata  et  al.  1994; 
Hansen  and  Hoppe  1997).  Figure  8  also  shows  the  reg¬ 
ular  downward  movement  of  wave  phase  fronts  above 
and  below  the  lamina,  where  the  response  approximates 
a  linear  wave  solution  of  the  form  (24). 

e.  Three-dimensional  gravity  wave  advection  of  a 

tracer 

The  features  that  produce  laminated  ozone  profiles 
tend  to  occur  in  long  narrow  streamers  of  meandering 
air,  which  lead  to  large  background  ozone  gradients  not 
only  vertically,  but  also  horizontally  (e.g.,  Waugh  et  al. 


1994;  Orsolini  et  al.  1997).  These  gradients  can  be  im¬ 
portant  when  attempting  to  isolate  wave-related  struc¬ 
tures  in  ozone  measured  within  such  environments  (e.g., 
Danielsen  et  al.  1991;  Teitelbaum  et  al.  1996;  Bac- 
meister  et  al.  1997;  Gibson- Wilde  et  al.  1997).  Again, 
the  approximate  relations  (24)  and  (28)  cannot  accu¬ 
rately  simulate  the  effects  of  gravity  waves  on  such 
three-dimensional  tracer  distributions. 

The  parcel-based  theory  in  appendix  A  assumed  no 
horizontal  gradients  in  the  background  tracer  profile  and 
cannot  be  directly  applied  here.  However,  the  analysis 
is  easily  extended  to  accommodate  full  three-dimen¬ 
sional  wave-induced  displacements  [^'(x,-,  y;,  z;,  t), 
r]' (x,.,  t),  £'(x;,  yt,  Zi,  0]  °f  a  collection  of  parcels, 

i  (i  =  1,  2,  3,  •  •  • )  at  equilibrium  positions  (x,,  yf,  z;) 
within  an  arbitrary  three-dimensional  tracer  distribution 
q(x,  y,  z). 

The  mixing  ratio  response  (A19)  generalizes  to 
q[xt  +  x'(xh  yt,  Zi,  t),  y,  +  tj'(x„  y;,  zt,  t), 

Zi  +  C(Xi,  yt,  Zi,  0]  =  q(x„  yt,  Zi),  (29) 

which,  as  discussed  in  section  3b,  holds  for  hydrostatic 
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and  nonhydrostatic  waves.  If  we  can  assume  that  back-  responding  extension  of  (A21)  yields  the  number  den- 
ground  isentropic  surfaces  are  horizontal,  then  a  cor-  sity  response 


n  [x,  +  x'(x „  y„  Zi,  t),  y,  +  tj'Ct.  y„  z„  t),  z,  +  £'(*,■>  y„  z„  f)]  =  n  (xn  yh  zt)  exp 


-£'(■*„  yt,  Zj,  t ) 

y h  ( Zi ) 


(30) 


for  a  background  number  density  distribution  n(x,  y,  z). 

As  discussed  in  section  3b,  (30)  is  accurate  for  hy¬ 
drostatic  gravity  waves  only.  In  environments  where 
large  narrow  horizontal  gradients  in  background  poten¬ 
tial  temperatures/densities  also  arise,  additional  com¬ 
pression  and  expansion  occurs  as  parcels  are  horizon¬ 
tally  advected  by  the  wave.  In  this  situation,  not  only 
must  the  exponential  in  (30)  be  modified  to  account  for 
this  extra  horizontal  compression/expansion,  but  ex¬ 
tended  wave  equations  may  be  needed  that  incorporate 
concomitant  modifications  to  the  wave  motion  (e.g.. 
Lamb  and  Shore  1992).  Since  background  isentropes 
rarely  tilt  from  the  horizontal  by  more  than  a  few  de¬ 
grees,  then  (30)  should  adequately  describe  gravity 
wave-induced  perturbations  of  a  three-dimensional  dis¬ 
tribution  of  tracer  number  densities  in  most  circum¬ 
stances. 

It  should  be  noted  that  the  parcel  relations  (29)  and 
(30)  are  Lagrangian  expressions.  A  typical  application, 
then,  would  be  in  supplementing  synoptic-scale  parcel 
advection  models  with  oscillatory  gravity  wave  advec- 
tion  effects  (e.g..  Pierce  et  al.  1994).  The  Eulerian  tracer 
perturbations  that  such  simulations  will  produce  (e.g., 
fluctuations  in  vertical  tracer  profiles)  will  be  quite  com¬ 
plicated  in  general. 


f.  Gravity  wave  perturbations  of  shorter-lived 
constituents 


The  parcel-based  approach  can  also  be  used  for  short¬ 
er-lived  constituents.  We  illustrate  this  by  using  it  to 
reproduce  some  well-known  results  for  ozone. 

While  the  lower-stratospheric  ozone  in  Figs.  7-8  acts 
as  a  tracer,  the  ozone  production/loss  term  in  (lb), 
R0  ,  becomes  nonnegligible  above  ~25-30  km.  Assum¬ 
ing  Chapman  photochemistry,  then 


R 


o3 


2  J o2no2 


2ko3Joflb3 

ko2nMn02  ’ 


(31) 


where  nx  is  the  number  density  of  constituent  x,  J x  is 
its  photolysis  rate,  and  kx  is  its  recombination  rate  (Bras- 
seur  and  Solomon  1984).  As  shown  in  appendix  C,  ap¬ 
plying  a  parcel-based  analysis  to  (2)  and  (31)  yields  a 
linearized  response  of  the  photochemical  term  RQ  to  an 
adiabatic  vertical  displacement  of  the  parcel  of  the  form 


=  - B(z2)q'0fz2 )  -  C  (z2)T'(z2),  (32) 

dt 

for  small  vertical  displacements  £,  where  q0j(z2)  = 
q0fz2 )  +  q'ofZi)  from  (A19)  and  T'(z2)  is  the  TMS 
temperature  perturbation  (A13).  Since  R0fz)  ~  0,  then, 
for  a  linear  q(z)  profile. 


dq0fz2) 

dt 


w(z2), 


(33) 


where  w(z2)  =  d£/dt  is  the  vertical  velocity  of  the  parcel 
at  z2  (see  Fig.  2a).  Thus  (32)  and  (33)  yield 


dq'ofzf) 

dt 


~B  {z2)q'(h{z2) 

-  C(z2)T'(z2).  (34) 


This  equation  was  first  derived  by  Lindzen  and  Goody 
(1965)  using  a  linearized  perturbation  expansion  of  (2) 
and  (31). 

We  now  substitute  the  wave  solution  (14).  If  l  is  small 
and  the  vertical  variation  of  background  quantities  in 
(34)  is  gradual  (as  assumed  in  appendix  C)  and  ap¬ 
proximately  linear  over  the  interval  zx  ±  l,  then  explicit 
parcel  transport  to  z2  can  be  removed  so  that  the  lin¬ 
earized  mixing  ratio  and  temperature  perturbations  in 
(34)  can  also  be  expressed  as  wave  solutions  of  the  form 
(14)  and  (16).  Substitution  of  those  solutions  into  (34) 
yields,  with  the  help  of  (21),  (A5),  and  (13),  the  peak 
amplitude  relation 


4oMi) 


dqofz^/dz  +  C(Zi)\ 
dT(z,)/dz  +  ra  M  ) 

x  f(ztl  (35) 


where  T„  is  the  dry-adiabatic  lapse  rate.  This  is  the  basic 
ozone-temperature  relation  derived  for  Kelvin  waves 
by  Randel  (1990).  It  illustrates  the  two  limiting  cases 
ofthe  response.  For  o>2  »  B2  and  ( Clop1  <C  [(dqQJdz)/ 
( dT/dz  +  r„)]2,  terms  on  the  right  of  (34)  are  small 
while  terms  on  the  left  are  large  and  dominate,  so  (35) 
approximates  the  tracer  solution  (20).  For  w2  <C  B2  and 
(C/co)2  »  [(dq0]/dz) /(dT/dz  +  Tu)]2,  photochemistry 
dominates  so  that  terms  on  the  right-hand  side  of  (34) 
are  much  larger  than  the  parcel  transport  terms  on  the 
left.  In  this  latter  case,  (35)  tends  toward  the  photo- 
chemistry-dominant  limit 
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Fig.  6.  Time-height  surface  rendering  of  the  perturbed  Na  number 
density  response  in  Fig.  5.  Both  plots  span  one  wave  period.  Panel 

(a)  shows  variations  over  the  full  height  range  of  the  Na  layer,  while 

(b)  focuses  on  the  SSL. 


(36) 

which  also  follows  on  setting  the  left-hand  side  of  (34) 
to  zero. 

This  parcel-based  analysis  assumes  adiabatic  parcel 
motion,  whereas  the  perturbation  (/'()  (z2)  also  produces 
a  diabatic  heating/cooling  perturbation  in  the  potential 
temperature  equation.  However,  Zhu  and  Holton  (1986) 
show  that  this  leads  to  very  minor  departures  from  adi- 
abaticity  for  hydrostatic  gravity  waves.  Thus  parcel- 
based  analyses  also  work  well  in  simulating  (hydro¬ 
static)  gravity  wave-induced  perturbations  of  constit¬ 
uents  with  short  chemical  lifetimes,  so  long  as  the  per¬ 
turbed  chemical  response  does  not  lead  to  significant 
diabatic  damping  of  the  wave;  otherwise,  a  linearized 


perturbation  analysis  should  be  pursued  (e.g.,  Zhu  and 
Holton  1986). 


4.  Summary  and  conclusions 

This  study  has  formalized  the  parcel  advection  meth¬ 
od,  applied  it  to  gravity  waves,  and  compared  its  find¬ 
ings  with  standard  results  from  linearized  perturbation 
expansions  of  the  governing  equations.  The  combined 
results  of  section  3  have  shown  that  the  two  approaches 
are  equivalent  in  most  circumstances.  However,  we  also 
identified  certain  situations  in  which  one  or  the  other 
method  gave  more  accurate  results. 

The  parcel-based  method  becomes  inaccurate  when 
simulating  perturbations  of  tracer  density  profiles  by 
nonhydrostatic  gravity  waves.  For  shorter-lived  con¬ 
stituents,  the  parcel  method  also  becomes  inaccurate  if 
the  wave-induced  photochemical  response  produces  a 
significant  diabatic  feedback  on  the  wave.  In  these  in¬ 
stances,  the  results  from  traditional  linearized  pertur¬ 
bation  analysis  should  be  used  instead. 

On  the  other  hand,  the  parcel-based  approach  is  the 
better  method  for  modeling  hydrostatic  gravity  wave 
perturbations  of  tracers  with  sharp  spatial  gradients  in 
their  background  distribution.  We  illustrated  this  by  us¬ 
ing  it  to  model  wave-induced  modulations  of  sporadic 
sodium  layers  and  ozone  laminae  (Figs.  5-8). 

While  we  have  restricted  attention  here  to  gravity 
waves,  the  parcel  methodology  described  here  can  also 
be  applied  to  other  types  of  waves,  such  as  free  and 
forced  Rossby  waves,  mixed  Rossby-gravity  waves, 
Kelvin  waves,  and  tides.  For  some  of  these  waves,  both 
meridional  and  vertical  wave  displacements  are  neces¬ 
sary  to  simulate  their  effects  on  tracers  (e.g.,  Randel 
1993;  Stanford  and  Ziemke  1993).  Parcel  advection 
models  have  already  been  used  extensively  to  simulate 
the  effects  of  stationary  Rossby  waves  on  stratospheric 
tracer  distributions. 
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(a)  cot  =  0° 


(b)  cot  =  60" 


(c)  cot  =  120" 


(d)  cot  =  180” 


(e)  cot  =  240° 


(f)  cot  =  300” 


Fig.  7.  As  in  Fig.  5  but  using  a  representation  of  the  ozone  layer  based  on  the  data  of  Gibson-Wilde  et  al.  (1997),  with  a  lamina  at  14 
km.  Thick  light  curves  show  the  unperturbed  profile.  Solid  curves  show  perturbed  profiles  due  to  a  gravity  wave  of  A,  =  2.3  km  and  l  = 
200  m. 


APPENDIX  A 


Fig.  8.  As  in  Fig.  6a  but  showing  the  ozone  number  density  re¬ 
sponse  from  Fig.  7. 


Perturbation  Formulas  Using  the  Parcel 
Methodology 

a.  Background  profiles 

The  ideal  gas  law  and  the  altimeter  equation  (11)  yield 
the  standard  background  density  relation 


T{Zl) 

P(z2)  =  p(z,)=—  exp 
T  (z2) 


ta 


Pfo)  exp 


~C 

Hp(z ,) 


(Al) 


where 


Hp(z i)  -  tffa)^l  +  (A2) 

is  the  density  scale  height.  On  invoking  (5)  and  (6)  at 
Zi  and  z2,  then,  with  the  aid  of  (11),  (Al)  and  (A2), 
background  potential  temperature  varies  with  height  as 
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7=r ,  T(z2)  (z2  -  Zi)(T  -  1) 


©  fe)  =  @  (Zi)=— —  exp 
T  (Zi) 


yH(zt) 


@(Zi)  exp 


£ 


HM’ 


(A3  a) 


(A3b) 


while  background  potential  density  varies  with  height 
as 


7S7  a  TV/  x P fe)  (fe-Zi) 
£>(z2)  =  D(Zi)— -  exp1 


P  (z, )  \ 


D(Zi)  exp 


tf„(z,)/’ 


(A4a) 

(A4b) 


where 


„  ,  x  dH(Zl) 

Hd(z ,)  =  tffo) -  H - - - 

y  dz 


8 


N2(zt) 


(A5) 


b.  Temperature  perturbations 

Figure  1  a  constructs  the  perturbation  of  potential  tem¬ 
perature  produced  by  the  adiabatic  vertical  displacement 
f  of  an  air  parcel  from  its  equilibrium  position  z, ,  as 
discussed  in  section  2b.  Since  ©  in  the  parcel  is  con¬ 
served,  the  displacement  produces  a  perturbation  from 
the  background  profile  at  height  z2  of 

®'(z2)  =  ©fo)  -  0(z2).  (A9) 

From  (A3b)  and  (A5),  this  implies  a  normalized  per¬ 
turbation: 


®'(z2)  _  (  -£  \  _  . 

©  fe)  expUDfei)j 

^  ~£  =  -A2(z.) 

Hoizi)  g 

=  _1  30  (z,). 

©  (Zj)  dz 


(AlOa) 


(AlOb) 


is  the  potential  density  scale  height,  g  is  gravitational 
acceleration,  and  N{Zl)  is  the  background  Brunt-Vaisala 
frequency  at  height  z,  ■ 

From  (8),  analogous  relations  for  the  background  po¬ 
tential  number  density  of  a  given  passive  trace  constit¬ 
uent  are 


where 


s  (z2-z i) 

v(z2)  =  v(Zl)—~  exp' 


p(Zi)  exp 


n  (Zi)  \  yH(zt) 
-£ 


HXzjy 


(A6a) 

(A6b) 


HXz,) 


H(Zl)  dn  (Z] )\ 

n(z  i)  dz  ) 


(A7) 


is  the  potential  tracer  density  scale  height.  The  simpli¬ 
fied  relation  (A6b)  is  accurate  when  height  variations 
of  the  n(zl )  profile  are  approximately  exponential  over 
the  interval  £.  For  small  £  (A6b)  is  a  good  approxi¬ 
mation  whenever  n(z)  varies  over  a  typical  vertical 
length  scale  L  »  £  However,  while  background  tem¬ 
peratures  in  (Al)  and  background  densities  in  (A4a) 
usually  obey  such  conditions,  n(z)  profiles  in  (A6a)  of¬ 
ten  do  not  (see  sections  3c  and  3d),  whereupon  (A6a) 
or  the  integrated  expression 


v(z2)  =  v(Zi)  exp 


(A8) 


must  be  used.  Note  that  the  reformulation  of  (A6a)  into 
a  fully  exponential  form  (A6b)  and  (A8)  is  similar  to 
the  approach  taken  by  Lindzen  (1981)  in  parameterizing 
the  variation  with  height  of  WKB  gravity  wave  equa¬ 
tions  [see  his  Eqs.  (11)  and  (12)]. 


where  (A  10b)  follows  from  retaining  the  first  two  terms 
of  the  MacLaurin  series  expansion  of  the  exponential 
in  (AlOa)  and  is  accurate  for  small  £  (i.e.,  \£/HD(Zl)\  <C 
1).  This  simplification  is  referred  to  hereafter  as  a  trun¬ 
cated  MacLaurin  series  (TMS)  expansion  and  from 
(AlOb)  is  clearly  accurate  when  the  background  profile 
@(z)  is  approximately  linear  between  z,  and  z2. 

The  absolute  temperature  response  is  depicted  in  Fig. 
lb.  Using  (5)  and  (11),  the  temperature  of  the  displaced 
parcel  becomes 

7-fa)  =  ffe,)  °*p(~^(~)1>)  (Alla) 

~r(z  i)-|F,  (Al  lb) 

'-'p 

where  Cp  is  the  mass  specific  heat  at  constant  pressure. 
Equation  (Allb)  follows  from  a  TMS  expansion  of 
(Alla)  and  reproduces  the  familiar  result  that  the  parcel 
temperature  varies  according  to  the  ( dry)  adiabatic  lapse 
rate  Fo  =  g/Cp  for  small  £  This  yields  a  temperature 
perturbation 

T\z2)  =  T(z2)  ~  T(za)  (A12) 

and  a  relative  perturbation,  according  to  (A5),  (Alla), 
and  (A12),  of 

Hz2)  _  Tfo)  (~£{y~  1)\  _  . 

T(z2)  T{z2)  expy  yH(Zl)  ) 


-  1 


~£  =  -A^(z,) 

Hd(Zi)  g 


(A13) 


using  a  TMS  expansion  of  the  exponential.  Thus  the 
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relative  perturbations  of  absolute  and  potential  temper¬ 
ature  are  equal  for  small 


c.  Density  perturbations 

Figure  lc  constructs  the  perturbation  of  potential  den¬ 
sity  produced  by  the  same  adiabatic  vertical  displace¬ 
ment  £  of  an  air  parcel  from  its  equilibrium  position  z, . 
Since  D  in  the  parcel  is  conserved,  then,  as  in  (A9), 

D'(z2)  =  D(z,)  -  D(Za).  (A14) 

From  (A4b)  and  (A5),  (A  14)  implies  a  normalized  per¬ 
turbation, 


D(z2) 


-  1 


i  N 2fa) 

HD(zt)  g  ’ 


(A15a) 

(A15b) 


where  (A15b)  follows  from  a  TMS  expansion  of  (A15a). 

Figure  Id  shows  the  corresponding  construction  for 
the  density  perturbation.  Using  (6)  and  (11),  the  density 
of  the  adiabatically  displaced  parcel  responds  according 
to 


Pfe) 


P(z  i)  exp 


yH{zx) 


(A  16) 


which,  by  analogy  to  the  adiabatic  lapse  rate  for  tem¬ 
perature  (A lib),  we  refer  to  here  as  the  adiabatic  rare¬ 
faction  curve.  This  yields  a  density  perturbation 


p'(z2)  =  p(z2)  -  p(z2), 

(A17) 

which,  with  the  aid  of  (Al),  (A2),  (A5), 
yields  to  a  relative  density  perturbation  of 

and  (A16), 

P'(z2)  _  p(z2 ) 
p(z2)  p(z2) 

‘  expUfe.))  ■ 1 

(A  18a) 

i  N2(Z,)  „ 

Hd{Z\)  g 

(A18b) 

using  a  TMS  expansion  of  (A18a). 


d.  Tracer  perturbations 

Figure  2a  constructs  the  perturbation  mixing  ratio  q' 
induced  by  the  adiabatic  vertical  displacement  £  of  a 
parcel  from  its  equilibrium  height  z,  in  an  atmosphere 
containing  a  tracer  with  a  background  vertical  profile 
q(z).  Since  q  in  the  parcel  is  conserved  according  to 
(3),  then 


q'(z2)  =  q(z2)  -  q  (z2) 
=  9(Zi)  “  q(z2) 


and 

q'(z2)  _  pUi)  _ 
q(z2)  q  (j.2) 

Since  p  =  MnM,  then  (3b)  and  (A16)  imply 


n(z2)  =  q  (z,  )nM{z.2) 


(A19) 


(A20) 


'  “p(isfc)'  <A2I) 

where  n(z)  is  the  background  number  density  profile  of 
the  tracer.  Thus,  like  the  total  density  p  in  (A16),  the 
tracer  number  density  n  responds  to  adiabatic  vertical 
displacements  according  to  its  adiabatic  rarefaction 
curve  (A21).  As  shown  in  Fig.  2b,  this  produces  a  per¬ 
turbation  of  tracer  number  density 

n'{z2 )  =  n(z2)  -  n(z2),  (A22) 

whose  relative  magnitude  is  given  by 

n'(z2) 
n  (z2) 


where  (A23c)  is  derived  from  (A23b)  using  (A  18a)  and 
(A20). 

An  equivalent  derivation  in  terms  of  potential  tracer 
number  density  is  shown  in  Fig.  2c.  Since  v  in  the  parcel 
is  conserved,  then 

v'  (z2)  =  v(z,  )  ~  v(z2),  (A24) 

which,  with  the  aid  of  (8),  (11),  and  (A23a),  yields 


v'{z,2) 

v{z2) 

_  viz,) 
v(z2) 

(A25a) 

«(m)  ~£  \  , 

n(,z2)  CXP\yW(Zi)/ 

(A25b) 

_  n’(z2) 

n  (z2) ' 

(A25c) 

If  £  is  small  and  v(z)  variations  are  of  much  longer 
vertical  scale  than  £  then  (A25c)  can  be  simplified  using 
a  TMS  expansion  of  (A6b),  along  with  (A7)  and  (A25a), 
to  yield 


n(z2) 

n{z2) 


-  1 


«fei)  (  ~j 

n(z2)  eXP\y//(z,) 

q(zi)  (  C 

q(z2 )  eXPUDfc) 


q'(z2)  p\z2)  (q'{z2)\(p'{z2) 


q(z2 )  pfe)  \q(z2)  \p(z2) 


(A23a) 

(A23b) 

(A23c) 
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Fig.  Bl.  Vertical  wavelength  variation  of  <zmax(m)  and  bmax(m )  (a), 
and  deviations  in  magnitude  (b)  and  phase  (c)  between  the  acoustic- 
gravity  wave  polarization  relation  (appendix  B)  and  the  corresponding 
parcel-based  approximation  (appendix  A)  relating  £  to  perturbations 
in  relative  temperature  (dashed  line)  and  relative  density  (dotted  line). 


v'(z2)  =  n'(z2)  ^  £ 

v(z2)  n  (z2)  HJz.,  ) 


(A26) 


The  same  expression  also  results  from  a  direct  TMS 
expansion  of  (A23a).  Note  the  similar  forms  of  (A26) 
and  (A15b). 


the  following  perturbation  equations  govern  a  zonally 
propagating  nonhydrostatic  gravity  wave  in  a  compress¬ 
ible  rotating  unsheared  inviscid  atmosphere  (e.g.,  Gos- 
sard  and  Hooke  1975): 

—  lcoCi  —  fv  =  —  ck&,  (B3) 

—  ccov  +  fu  =  0,  (B4) 

—  lcow  +  rg  =  —  (lM  —  (B5) 


—  Lcor  +  iku  +  (iM  —  H  pl)w  =  0,  (B7) 

where  cs  is  background  sound  speed,  (n,  v,  w)  is  the 
peak  velocity  amplitude,  r  =  p/p  is  the  relative  density 
amplitude,  and  (\>  =  pl~p ,  in  which  p  is  the  pressure 
amplitude.  The  full  nonhydrostatic  dispersion  relation 
for  gravity  waves  [e.g.,  (23-7)  of  Gossard  and  Hooke 
(1975)]  follows  from  the  equation  set  (B3)-(B7);  that  is. 


m 2  =  k 2 


1  uE 

~Aw-p  +  y 


(B8) 


Polarization  relations  follow  from  (B3)-(B7).  Elimi¬ 
nating  between  (B5)  and  (B6)  and  noting  that  w  = 
—  lw£.  then 


[1  +  a  (m,  co)  +  ib(m,  cj)]£. 


(B9) 


a(m,  oj) 


(1  -  co2/N2) 


(B10) 


b(m,  co) 


(1  -  co2/N2) 


m 

+  r- 


(B 11) 


where  T  =  —1/(2 Hp)  +  g/c2  is  the  anisothermal  form 
ofEckart’s  coefficient.  Since  c2  =  ygH,  then  eliminating 
r  between  (B5)  and  (B6)  yields 


<f>  =  N2yH[a(m,  co)  +  cb(m,  co)]l,  (B12) 


where  y  is  the  ratio  of  specific  heats.  This  replicates 
(23-6)  of  Gossard  and  Hooke  (1975)  and  can  be  reex¬ 
pressed  in  terms  of  the  relative  pressure  perturbation 


p  _  d> 
p  ~  RT 


y [a  (in,  co)  +  ib(m,  co)]£. 


(B13) 


APPENDIX  B 

Acoustic-Gravity  Wave  Polarization  Relations 

On  splitting  atmospheric  field  variables  into  the  sum 
of  a  mean  and  perturbation  term  (X  =  X  +  X'),  isolating 
and  linearizing  perturbation  terms,  and  assuming  wave 
solutions  of  the  form 

X'  =  X  exp[i(kx  +  Mz  —  Hf)]  (Bl) 

=  X  exp[i(kx  +  mz  —  Ht)  +  z!2Hp],  (B2) 


A  linearized  perturbation  expansion  of  the  ideal  gas  law 
yields 


f 

T 


P 

P 


(B14) 


which,  on  inserting  (B9)  and  (B13),  yields 


f 

T 


{ 1  —  (y  —  1  )[a(m,  co)  +  ib (m,  co)]}£. 

(B15) 


15  December  1998 


ECKERMANN  ET  AL. 


3537 


Equations  (B3)-(B7)  also  yield  the  polarization  relation 
-O2  -  f2) 

w  =  - (m  +  lT)u.  (B16) 

k(N2  -  w2) 


phase  from  the  full  acoustic-gravity  wave  polarizations 
relations  (B9)  and  (B 15). 


For  an  isothermal  ( H  =  Hp)  nonrotating  (/  =  0)  at¬ 
mosphere,  it  can  be  shown  that  the  dispersion  and  po¬ 
larization  relations  (B8)-(B16)  are  the  same  as  Eqs. 
(21)— (24)  and  (26)  of  Hines  (1960). 

On  linearizing  the  tracer  continuity  equation  (lb),  as¬ 
suming  a  height-varying  mean  tracer  profile  n{z),  and 
substituting  the  wave  solution  (Bl),  we  obtain 

dn  _  _  - 

—  LU>n  H - w  +  n  V  •  U  =  0. 

dz 

From  (B16),  it  follows  that 
V  •  U  =  iku  +  lMw, 

f  1  N2  ,  1 

=  i - [a(m,  co)  +  Lb(m,  «)]  \w, 

\yn  8 


APPENDIX  C 

Photochemical  Reponse  of  Ozone  Mixing  Ratio  in 
a  Vertically  Displaced  Air  Parcel 

We  start  with  the  parcel-based  ozone  rate  equation 
(Lindzen  and  Goody  1965) 


dq0i 

(B17) 

dt 

where 

(B18) 

ko2 

(B19) 

2ek0iJ0}qd} 


300  K) 

X  exp [  —  (2060  K)/T]  =  G(T ) 


(Cl) 


(C2) 


whereupon  (B17)  yields 

n  f  1  1  dn  N2  ]  . 

-  -  \  —  +  - [a(m,  co)  +  cb(m,  co)]  \  l- 

n  I  yH  n  dz  g 

(B20) 

Thus,  on  comparing  (B9),  (B13),  (B15),  and  (B19) 
with  the  simplified  parcel-based  derivations  in  appendix 
A  (see  also  sections  2  and  3),  we  see  that  the  normalized 
terms  I a(m,  w)l  and  I b(m,  o>)l  must  be  <Cl  for  the  two 
sets  of  equations  to  be  equal  to  a  good  approximation. 
We  note  from  (B10)  and  (Bll)  that  the  frequency  de¬ 
pendence  of  a{m,  oj)  and  b(m,  co)  is  weak,  so  for  sim¬ 
plicity  we  consider  the  limiting  maximum  values 


-r 


+  r- 


(B21) 


g  I  m  \ 

b„rJm)  =  ^  2^rJ  (B22) 

c2\m2  +  1  -J 

that  arise  for  ( co/N )2  — »  0. 

Since  these  maximum  values  depend  only  on  the  ver¬ 
tical  wavenumber  m.  Fig.  Bla  plots  amdJm)  and  bmilJm) 
versus  A.  =  27r/lml  for  an  isothermal  atmosphere  of  T 
=  230  K.  Figures  Bib  and  Blc  show  the  corresponding 
maximum  difference  in  magnitude  and  phase  between 
the  polarization  relations  (B9)  and  (B15)  and_the  cor¬ 
responding  parcel  approximations  p/p  =  —TIT  =  ( N 2 1 
g)l  given  by  (18).  Other  plots  (not  shown)  revealed  that 
these  differences  increase  (decrease)  with  decreasing 
(increasing)  T,  and  that  the  equality  in  magnitude  of  the 
density  and  temperature  departures  in  Fig.  Bib  occurs 
only  when  the  atmosphere  is  isothermal.  For  vertical 
wavelengths  A_  S  20  km,  we  see  from  Fig.  Bl  that  the 
approximate  parcel-based  formulas  of  appendix  A,  ap¬ 
proximated  by  the  wave  polarization  relations  (14)  and 
(18),  produce  only  small  deviations  in  amplitude  and 


(Demore  et  al.  1994).  Equation  (Cl)  differs  from  (2) 
and  (31)  only  through  insertion  of  an  extra  factor,  e, 
that  scales  the  ozone  destruction  rate  due  to  recombi¬ 
nation  with  monatomic  oxygen,  kl2,  to  a  net  rate  that 
resembles  the  one  obtained  when  reactions  with  NO,, 
HOv,  and  CIO,  species  are  included  (Zhu  and  Holton 
1986). 

We  use  (Cl)  and  (C2)  to  describe  the  time  evolution 
of  ozone  in  a  parcel  displaced  adiabatically  by  a  small 
vertical  distance  Since  (A  18a)  holds  for  «0,(z2)/ 
nQl(z2 )  -  1,  then 


dq0,(z2) 

dt 


27o2p  0,(^2) 


270, 

«o2fe) 


G(T(z2))p63(z2)  exp (~C/Hd), 

(C3) 


where  q0l(z2)  =  qQi{z2)  +  qo3(z2)  from  (A19).  We  sim¬ 
plify  (C3)  by  first  noting  that  the  background  produc¬ 
tion/loss  rate 


3fe)  2  J Q}q  o,  (z2 )  27, 


G(T(z2))q23(z2) 


0. 


(C4) 


We  then  use  a  TMS  expansion  of  exp(— £7 HD)  and  a 
truncated  Taylor  series  expansion 


G(T  +  T')  =  G(T)  +  T'{^j  =  G  +  G',  (C5) 

where  T'(z2)  =  T(z2)  —  T(z2).  Finally,  on  assuming 
constant  qQd  (no  diffusive  separation),  constant  7,,  and 
retaining  only  linear  terms  in  the  perturbation  expansion 
of  variables  in  (C3),  we  obtain,  with  the  aid  of  (A13), 
(C4),  and  (C5), 
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dt 


=  ~B  fe)<?o:ife)  -  C(z2)T'(z2),  (C6) 


where 


fife)  = 


_  4J0,gfe)G(r(z2)) 


«o2fe) 


C(z2)  =  — = 


B(z2)q(z2) 


2T*(z2) 


T(z2 ) 


/  cf(lnG)\ 

U1/r>/r, 


(C7a) 

.  (C7b) 


Equation  (C6)  then  yields  the  well-known  ozone  per¬ 
turbation  result  of  Lindzen  and  Goody  (1965)  (see  sec¬ 
tion  3f),  with  the  expressions  (C7a)  and  (C7b)  essen¬ 
tially  the  same  as  those  derived  by  Zhu  and  Holton 
(1986).1  More  complete  expressions  for  RQ  with  the 
NO,.,  HOv,  and  CIO,  rate  terms  included  explicitly  also 
leads  to  a  linearized  ozone  deviation  equation  of  the 
same  form  as  (C6),  although  the  formulas  for  the  co¬ 
efficients  B  and  C  are  more  complicated  in  this  case 
(see,  e.g.,  Stolarski  and  Douglass  1985). 
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